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We compute the beta functions for the three gauge couplings of the Standard 
Model in the minimal subtraction scheme to three loops. We take into account 
contributions from all sectors of the Standard Model. The calculation is performed 
using both Lorenz gauge in the unbroken phase of the Standard Model and back- 
ground field gauge in the spontaneously broken phase. Furthermore, we describe in 
detail the treatment of 75 and present the automated setup which we use for the 
calculation of the Feynman diagrams. It starts with the generation of the Feynman 
rules and leads to the bare result for the Green's function of a given process. 
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I. INTRODUCTION 



Renormalization group functions are fundamental quantities of each quantum field theory. 
They provide insights in the energy dependence of cross sections, hints to phase transitions 
and can provide evidence to the energy range in which a particular theory is valid. The 
renormalization group functions of the gauge couplings in the Standard Model (SM) are of 
particular importance in the context of Grand Unified Theories allowing the extrapolation 
of low-energy precision data to high energies, not accessible by collider experiments. 

Important milestones for the calculation of the gauge coupling beta functions in the 
Standard Model are the following computations: 

• The one-loop beta functions in gauge theories along with the discovery of asymptotic 
freedom have been presented in Refs. [j], 0] • 

• The corresponding two-loop corrections 

— in gauge theories without fermions [H, H[ , 

— in gauge theories with fermions neglecting Yukawa couplings [5-7], 

— with corrections involving Yukawa couplings jsf, 

are also available. 



The two-loop gauge coupling beta functions in an arbitrary quantum field theory have 
been considered in Ref. 0, 10 . 



The contribution of the scalar self-interaction at three-loop order has been computed 
in 



11, 12 



The gauge coupling beta function in quantum chromodynamics (QCD) to three loops 
is known from Ref. 13|, [14 . 



The three-loop corrections to the gauge coupling beta function involving two strong 
and one top quark Yukawa coupling have been computed in Ref. 
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The three-loop corrections for ageneral quantum field theory based on a single gauge 
group have been computed in [16 . 



The four-loop corrections in QCD are known from Refs. [17J, [18 . 



Two-loop corrections to the renormalization group functions for the Yukawa and Higgs boson 
self-couplings in the Standard Model are also known (loL 19-23]. Recently the dominant 
three-loop corrections to the renormalization group functions of the top quark Yukawa and 
the Higgs boson self-coupling have been computed in Ref. [24| . In this calculation the gauge 
couplings and all Yukawa coupling except the one of the top quark are set to zero. 

In this paper we present details to the three-loop calculation of the gauge coupling renor- 
malization constants and the corresponding beta functions in the SM taking into account 



all sectors. The results have already been presented in Ref. [25 



The remainder of the paper is organized as follows: In the next Section we introduce our 
notation and describe in detail how we proceed to obtain the beta functions of the gauge 
couplings. In particular, we describe the calculation in Lorenz gauge and background field 
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gauge (BFG), discuss our setup for an automated calculation, and explain our treatment of 
7 5 . The analytical results for the beta functions are presented in Section IHIl In contrast 
to Ref. [25| we show the results including all Yukawa couplings. Section [IV] is devoted to a 
description of the checks which have been performed to verify our result. A discussion of the 
numerical impact of the newly obtained corrections is given in Section |V] We conclude in 
Section IVIl Explicit results for the renormalization constants are relegated to Appendix |A] 
and Appendix [Bl In Appendix O we present three- loop results for the beta functions of the 
QED coupling constant and the weak mixing angle. Furthermore, we present in Appendix [D] 
translation rules which are useful in order to compare parts of our findings with the results 
of Ref. 0. 



II. THE CALCULATION OF THE BETA FUNCTIONS 

In this paper we present the beta functions for the three gauge couplings of the SM 
up to three loops in the modified minimal subtraction (MS) renormalization scheme. In 
the corresponding calculation we took into account contributions involving the three gauge 
couplings of the SM, the top, the bottom, and the tau Yukawa couplings and the Higgs 
self-coupling. Let us mention that we were able to derive the beta functions for a general 
SM Yukawa sector from the calculation involving the aforementioned seven couplings. We 
postpone the discussion of this issue to the next Section. In the following, we give details 
on the computation at the three-loop order. 

We define the beta functions as 

^-ftifeM, (i) 

where e = (4 — d)/2 is the regulator of Dimensional Regularization with d being the space- 
time dimension used for the evaluation of the momentum integrals. The dependence of the 
couplings ai on the renormalization scale is suppressed in the above equation. 

The three gauge couplings a±, and 03 used in this paper are related to the quantities 
usually used in the SM by 

5 CtQED 



Q 2 



3 cos 2 6 W 

«QED 



sin 9 W 

a 3 = a s , (2) 

where «qed is the fine structure constant, 6 W the weak mixing angle and a s the strong 
coupling. We adopt the SU(5) normalization which leads to the factor 5/3 in the equation 
for Gi\. 

The Yukawa couplings are defined as 

a *=2^e w M* w Wlth x = M ' r ' (3) 

where m x and Mw are the fermion and W boson mass, respectively. 
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We denote the Higgs boson self-coupling by A, where the Lagrange density contains the 
following term 



£sm = •••" M) (H^H) 2 + ... , (4) 



describing the quartic Higgs boson self-interaction. 

The beta functions are obtained by calculating the renormalization constants relating 
bare and renormalized couplings via 



where i = 1, 2 or 3. We furthermore set a 4 = a t , a 5 = a e = a T and a 7 = A. 

The first term in the first factor of Eq. ([6]) originates from the term /i 2<E in Eq. (jSJ) and 
vanishes in four space-time dimensions. The second term in the first factor contains the beta 
functions of the remaining six couplings of the SM. Note that (for the gauge couplings) the 
one-loop term of Z a . only contains a iy whereas at two loops all couplings are present except 
A. The latter appears for the first time at three-loop level. As a consequence, it is necessary 
to know (3j for j — 4, 5, 6 to one- loop order and only the e-dependent term for (3?, namely 
/?7 = —60:7/77. From the second term in the first factor and the second factor of Eq. ([6]) one 
can read off that three-loop corrections to Z ai are required for the computation of Pi to the 
same loop order. 

We have followed two distinct paths to obtain the results for the three-loop renormaliza- 
tion constants, which we discuss in the following Subsections, where we mention features 
and differences. 



The first method used for the calculation of the renormalization constants is based on 
Feynman rules derived for the SM in the unbroken phase in a general Lorenz gauge with 
three independent gauge parameters corresponding to the three simple gauge groups. All 
building blocks of our calculation are evaluated for general gauge parameters in order to 
have a strong check of the final results for the beta functions which have to be gauge 
parameter independent. It is possible to use the unbroken phase of the SM since the beta 
function in the MS scheme is independent of all mass parameters and thus the spontaneous 
symmetry breaking does not affect the final result. Note that this choice is advantageous 
for the calculation because in the unbroken phase much less different types of vertices have 
to be considered as compared to the phase in which the spontaneous symmetry breaking is 
present. 

In principle each vertex containing the coupling g { = ^Aixa>i can be used in order to 
determine the corresponding renormalization constant via 




(5) 



Taking into account that a} are does not depend on /1, Eqs. ([I]) and fl5]) lead to 




(6) 



A. Lorenz gauge in the unbroken phase of the SM 



Z. 




(7) 
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where Z vvtx stands for the renormalization constant of the vertex and Z^^t for the wave 
function renormalization constant; k runs over all external particles. 

We have computed Z a3 using both the ghost-gluon and the three-gluon vertex. Z a2 has 
been evaluated with the help of the ghost- W3, the W1W2W3 and the <fi + (fi~W3 vertex where 
± is the charged component of the Higgs doublet corresponding to the charged Goldstone 
boson in the broken phase and Wi, W2 and W3 are the components of the W boson. As to 
Z ai , a Ward identity guarantees that there is a cancellation between the vertex and some 
of the wave function renormalization constants yielding 

^ai = (8) 

where Zb is the wave function renormalization constant for the gauge boson of the U(l) 
subgroup of the SM in the unbroken phase. 

In Fig. [I] we show several one-, two- and three-loop sample diagrams contributing to the 
considered two- and three-point functions. 

We have not used vertices involving fermions as external particles as they may lead to 
problems in connection with the treatment of 75 in d 7^ 4 dimensions. The vertices selected 
by us are safe in this respect. A detailed discussion of our prescription for 75 is given below. 

In order to compute the individual renormalization constants entering Eq. (jjTj) we proceed 



as outlined, e.g., in Ref. [15]. The underlying formula can be written in the form 

Z r = 1 - K t (Z r T) , (9) 

where T represents the two- or three-point function corresponding to the renormalization 
constant Zr and the operator K e extracts the pole part of its argument. From the structure 
of Eq. (jHJ) it is clear that Zr is computed order-by-order in perturbation theory in a recursive 
way. It is understood that the bare parameters entering F on the right-hand side are replaced 
by the renormalized ones before applying K e . The corresponding counterterms are only 
needed to lower loop orders than the one which is requested for T. In our approach the three- 
loop calculation of Z ai requires — besides the result for Z ai to two loops — the one-loop 
renormalization constants for the other two gauge and the Yukawa couplings. Furthermore 
we have to renormalize the gauge parameters; the corresponding renormalization constants 
are given by the wave function renormalization constants of the corresponding gauge bosons 
which we anyway have to evaluate in the course of our calculation. 



B. Background field gauge in the spontaneously broken phase 

The second method that we used in order to get an independent result for the renormal- 



ization constants of the gauge couplings is a calculation in the BFG [26|, |27|. The basic idea 
of the BFG is the splitting of all gauge fields in a "quantum" and a "classical" part where 
in practical calculations the latter only occurs as external particle. 

The BFG has the advantage that Ward identities guarantee that renormalization con- 
stants for gauge couplings can be obtained from the exclusive knowledge of the corresponding 
wave function renormalization constant Q Thus we have the following formula 



1 In Lorenz gauge, this only works for £7(1) gauge groups, cf. Subsection III Al 
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FIG. 1. Sample Feynman diagrams contributing to the Green's functions which have been used 
for the calculation of the renormalization constants of the gauge couplings. Solid, dashed, dotted, 
curly and wavy lines correspond to fermions, Higgs bosons, ghosts, gluons and electroweak gauge 
bosons, respectively. 
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Zo, = ^— , (10) 

where A denotes the gauge boson corresponding to the gauge coupling a^. 

In contrast to the calculation using Lorenz gauge, we performed the calculation in the 
BFG in the spontaneously broken phase of the SM. As discussed in the last Subsection, such 
a calculation is more involved than a calculation in the unbroken phase since more vertices 
are present. On the other hand it constitutes an additional check of our result, that allows 
us not only to compare the BFG and the Lorenz gauge but also to switch from the broken 
to the unbroken phase of the SM . 

Since the calculation has been performed in the broken phase we have computed the 
transverse part of the two-point functions of the (background) photon, Z boson, photon- Z 
mixing, W boson and gluon which we denote by Il 7 , Hz, II 7 z, H-w, n s , respectively. Sample 
Feynman diagrams up to three loops can be found in the first two lines of Fig. [TJ 

H\y and U g can be used in analogy to Subsection III Al in order to obtain the corresponding 
renormalization constants which leads in combination with Eq. fllOp to the renormalization 
constants for a 2 and a s . We found complete agreement with the calculation performed in 
Lorenz gauge. 

As far as the self energies involving photon and Z boson are concerned, we consider at the 
bare level appropriate linear combinations in order to obtain the B and W boson self-energy 
contributions. To be precise, we have 

U B = cos 2 ^ re IT 7 + 2 cos ^ re sin 9^ e U lZ + sin 2 6^ re U z , 

LV = sin 2 ^ re n 7 - 2 cos^ re sin 6^ Te Tl lZ + cos 2 ^TIT Z . (11) 

The second linear combination can immediately be compared with the bare result obtained 
from the charged W boson self energy and complete agreement up to the three-loop order has 
been found. This constitutes a strong consistency check on the implementation of the BFG 
Feynman rules. LL3 is used together with Eq. (TTOT) in order to obtain the renormalization 
constant for a±. Again, complete agreement with the calculation described in the previous 
Subsection has been found. 

In our BFG calculation we want to adopt Landau gauge in order to avoid the renor- 
malization of the gauge parameters £j. However, it is not possible to choose Landau gauge 
from the very beginning since some Feynman rules for vertices involving a background gauge 
boson contain terms proportional to where & = corresponds to Landau gauge. To 
circumvent this problem we evaluate the bare integrals for arbitrary gauge parameters. In 
the final result all inverse powers of £j cancel and thus the limit £j = can be taken at the 
bare level. 

C. Automated Calculation 

Higher order calculations in the SM taking into account all contributions are quite in- 
volved. Apart from the complicated loop integrals, there are many different Feynman rules 
and plenty of Feynman diagrams which have to be considered. In our calculation we have 
used a setup which to a large extend avoids manual interventions in order to keep the 
error-proneness to a minimum. 
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FIG. 2. Overview of our automated setup. Calling up the programs in the uppermost line de- 
termines and evaluates a given process in a given model. The vertical workflow leads to the 
implementation of a new model in the setup. The programs are discussed in more detail in the 
text. 



As far as the loop integrals are concerned we exploit the fact that the beta function in 
the MS scheme is independent of the external momenta and the particle masses. Thus, we 
can choose a convenient kinematical configuration which leads to simple loop integrals as 
long as the infra-red structure is not modified. In our case we set all particle masses to 
zero and only keep one external momentum different from zero. We have checked that no 
infra-red divergences are introduced as we will discuss in detail in Section IIVI In this way 
all loop-integrals are mapped to massless two-point functions that up to three loops can be 



computed with the help of the package MINCER 28 



As core of our setup we use a well-tested chain of programs that work hand-in- 



hand: QGRAF [29[ generates all contributing Feynman diagrams. The output is passed via 



q2e 30, 31], which transforms Feynman diagrams into Feynman amplitudes, to exp [30|, 31 1 
that generates FORM j32[ code. The latter is processed by MINCER [28] and/or MATAD [33 j 
that compute the Feynman integrals and output the e expansion of the result. The paral- 
lelization of the latter part is straightforward as the evaluation of each Feynman diagram 
corresponds to an independent calculation. We have also parallelized the part performed by 
q2e and exp which is essential for our calculation since it may happen that a few times 10 5 
diagrams contribute at three-loop level to a single Green's function. The described workflow 
is illustrated on the top of Fig. [2J 

In order to perform the calculation described in this paper we have extended the above 
setup by the vertical program chain in Fig. [21 The core of the new part is the program 
FeynArtsToQ2E which translates FeynArts [34( model files into model files processable by 
QGRAF and q2e. In this way we can exploit the well-tested input files of FeynArts in our 
effective and flexible setup based on QGRAF, q2e, exp and MINCER. This avoids the coding of 
the Feynman rules by hand which for the SM would require a dedicated debugging process. 

For the part of our calculation based on the BFG we have used the FeynArts model files 
which come together with version 3.5. However, for Lorenz gauge in the unbroken phase 
there is no publicly available FeynArts model. For this reason we have used the package 
FeynRules [35[ in order to generate such a file which is also indicated in Fig. [2j 

Let us mention that FeynArtsToQ2E is not restricted to the SM but can process all model 
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files available for FeynArts. 

D. Treatment of 75 

An important issue in multi-loop calculations is the definition of 75 away from d — 4 
dimensions. A first possibility is the naive regularization that requires that 75 anti-commutes 
with all other 7-matrices. This approach has the advantage that its implementation is very 
simple. However, it can lead to wrong results, especially for Feynman diagrams involving 
several fermion loops. For example, the naive regularization of 75 leads to the problematic 
result (see, e.g., Ref. (36j |) 

tr( 7 "7<W75) = (d ^ 4) . (12) 

The limit of this expression for d — > 4 does not agree with its value in the physical case, 
when the regularization is turned off 

ti{YYYi a i5) = -&e» vp ° (d = 4) . (13) 

Here the totally anti-symmetric Levi-Civita tensor is defined by £ 0123 = 1. 

It is therefore reassuring that one can show explicitly that in the computation via the 
ghost-ghost-gauge boson vertex^] all contributions stemming from this kind of traces vanish. 
To prove this, we notice that this kind of traces can only lead to non- vanishing contributions 
if there are at least two of them in a diagram. Only in this case the e-tensors originating 
from Eq. (fT3"|) can be contracted, providing Lorentz structures that may contribute to the 
renormalization constants. We observe that the fermion loops can only yield problematic 
non-vanishing contributions if at least three lines are attached to them. Otherwise, there are 
too few external momenta and too few open Lorentz indices available. A general three-loop 
diagram with at least two closed fermion loops has the following form 




Here the solid lines represent fermions while the double lines can be either scalar bosons or 
gauge bosons. 

One can easily show that one cannot attach ghosts as external particles in the above dia- 
gram, since there are no vertices involving ghosts and fermions. Thus, the ghost self-energy 
and the ghost-ghost-gauge boson vertex do not contain diagrams of this type. Therefore, we 
need to consider only diagrams with two external gauge bosons for the following discussion. 
So the diagrams which still have to be discussed have the structure 




2 We restrict the discussion in this Subsection to the ghost-ghost-gauge boson vertex. For all other vertices 
without external fermions the reasoning is in complete analogy. 
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There are three ways to replace the double lines: 




The first diagram type cannot yield problematic contributions for the same reason as men- 
tioned above. The second diagram type vanishes as the fermion traces involve exactly five 
7-matrices, not counting 75 matrices. (Remember that we deal with diagrams in which all 
the propagators are massless.) Finally, the third diagram type involves fermion loops with 
three external gauge bosons. Such diagrams can indeed contain contributions originating 
from traces of 75 and an even number of other 7-matrices. However, the sum over all pos- 
sible fermion species that can circulate in the loops, including also the diagrams in which 
the fermions circle in opposite directions, vanishes, j This is of course a consequence of the 
cancellation of the Adler-Bell-Jackiw anomaly [38|, 3^] within the SM, as required by gauge 
invariance. Therefore, we are allowed to calculate the Feynman diagrams contributing to Z a . 
using a naive regularization prescription for 75, in which the diagrams containing triangle 
anomalies are set to zero from the very beginning, according to Eq. (IT21) . 

As an additional check of the calculation we implemented also a "semi-naive" regulariza- 
tion prescription for 75. Explicitly, we evaluate the expression tr^^^^^s) by applying 
the formal replacement 

tr(7^7V7 ff 7 5 ) = -^ VfXT + 0(e) . (14) 
The tensor e^p" has some similarities with the four-dimensional Levi-Civita tensor: {%) it is 
completely antisymmetric in all indices; {%%) when contracted with a second one of its kind 
one obtains the following result 



~livpcr ~ 



hi 

9^9, 



9n'9j 



(15) 



where the square brackets denote complete anti-symmetrization. When taking the limit 
d — > 4, e^ v P a converts into the four-dimensional Levi-Civita tensor and Eqs. (j!4p and ffl5|) 
ensure that it provides the correct four-dimensional result. 

At this point a comment on Eqs. (|14p is in order. It is straightforward to see that the 
combination of this equation and the cyclic property of traces leads to an ambiguity of order 
(9(e). Therefore, we made sure that the terms that need to be treated in this way generate at 
most simple poles in e and the above procedure can be applied directly without introducing 
additional finite counterterms. 

Let us stress again that we find the same result for the renormalization constants Z a . 
both from the ghost-ghost-gauge boson vertex and by using other vertices and both by 
applying the "naive" as well as the "semi-naive" scheme. These findings strongly support 
the above reasoning. 



E. Comparison of the methods 

This Subsection is devoted to a brief comparison of the calculation via the Lorenz gauge 
and the one involving the BFG. As has been mentioned before, in the BFG it is sufficient 



3 The proof of this statement can be based only on considerations about group theoretic invariants. For 
details see Chapter 20 of Ref. 37 1 . 
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Lorenz gauge BFG 



# loops 


1 


2 


3 


4 


# loops 


1 


2 


3 


4 


BB 


14 


410 


45 926 


7 111 021 




13 


416 


61 968 


13 683 693 


W 3 W 3 


17 


502 


55 063 


8438 172 


P, rypj 


13 


604 


100 952 


23 640 897 


99 


9 


188 


17611 


2 455 714 


Z Z 


20 


1064 


183 465 


44 049 196 


c g c g 


1 


12 


521 


46 390 


TT7-I- RTT7— P, 


18 


1438 


252 162 


42 423 978 


CW3CW3 


2 


42 


2 480 


251 200 


R B 

9 9 


10 


186 


17494 


2 775 946 


(f) + (/) 


10 


429 


46 418 


6 918 256 












BBB 


34 


2172 


358 716 


73 709 886 












w x w 2 w z 


34 


2 216 


382 767 


79 674 008 












999 


21 


946 


118 086 


20 216 024 












c g c g g 


2 


66 


4 240 


460 389 












c Wl c W2 W 3 


2 


107 


10 577 


1517 631 














24 


2 353 


387338 


77292 771 













TABLE I. The number of Feynman diagrams contributing to the Green's functions evaluated in 
this work. Left table: two- and three-point functions computed in Lorenz gauge; right table: 
two-point functions computed in BFG. The superscript "B" denotes background fields. The first 
column indicates the external legs of the Green's function, the other columns show the number of 
diagrams at the individual loop orders. Note that the BBB vertex is computed in order to have a 
cross check as we will explain in Section IIVI 

to consider only the gauge boson propagators. This is advantageous as the use of Lorenz 
gauge also requires the evaluation of three- (or four-) point functions and in most cases 
it also demands the consideration of additional two-point functions apart from the gauge 
boson ones. The disadvantages of the BFG are the increased number of vertices and the 
more involved structure of the vertices containing a background field. 

In Tab.HJwe list the number of diagrams for each Green's function contributing to the one-, 
two- and three-loop order. The number of diagrams computed in this work is obtained from 
the sum of the numbers in these columns. For comparison we also provide the corresponding 
number of diagrams which contribute to the four-loop order. 

It is tempting to compare the number of contributing Feynman diagrams in Lorenz gauge 
and in BFG which is, however, not straightforward since we use the former in the broken 
and the latter in the unbroken phase. Nevertheless, one observes that in the case of 03 the 
number of diagrams entering the BFG calculation is roughly the same as in case the gluon- 
ghost vertex is used in Lorenz gauge, even up to four-loop order. All other vertices lead 
to significantly more diagrams. In the case of P2 there are about three to four times more 
diagrams to be considered in the BFG as compared to Lorenz gauge. Whereas at three- 
loop order the difference between approximately 70 000 and 250 000 diagrams is probably 
not substantial it is striking at four-loop order where the number of Feynman diagrams 
goes from about 10 000 000 (VF3VF3 , cw 3 cw 3 and cw 1 cw 2 W 3 Green's function) to 42 000 000 
(W +B W~ B Green's function) when switching from Lorenz gauge to BFG. Thus, starting 
from four loops it is probably less attractive to use the BFG. 

Let us add that the precise number of Feynman diagrams depends on the detailed setup 
as, e.g., on the implementation of the four-particle vertices. Because of the colour structure 
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we split in our calculation the four-gluon vertex into two cubic vertices by introducing non- 
propagating auxiliary particles, whereas all other four-particle vertices are left untouched. 

Let us finally mention that the CPU time for the evaluation of an individual diagram 
ranges from less than a second to few minutes. For general gauge parameters it may take 
up to the order of an hour. Thus the use of about 100 cores leads to a wall-clock time which 
ranges from a few hours for a calculation in Feynman gauge up to about one day for general 
gauge parameters. For the preparation of the FORM files using QGRAF, q2e and exp also a few 
hours of CPU time are needed which is because of the large amount of Feynman diagrams. 
The use of about 100 cores leads to a wall-clock time of a few minutes. 



III. ANALYTICAL RESULTS 

In this Section we present the analytical results for the beta functions. As mentioned 
before, we are able to present the results involving all contributions of the SM Yukawa sector. 
The SM Yukawa interactions are described by (see, e.g., Chapter 11 of Ref. jioj ) 

Yukawa = -Qf Yg eH*uf - QfYgHdf - LfYjjHlf + h.c. , (16) 

where Y U,D,L are complex 3x3 matrices, i,j are generation labels, H denotes the Higgs 
field and e is the 2x2 antisymmetric tensor. Q L ,L L are the left-handed quark and lepton 
doublets, and u R , d R , l R are the right-handed up- and down-type quark and lepton singlets, 
respectively. The physical mass-eigenstates are obtained by diagonalizing Y U,D,L by six 
unitary matrices V^ ,L as follows 

YL g = VlYfvl\ f = U,D,L. (17) 

As a result the charged-current W ± couples to the physical quark states with couplings 
parametrized by the Cabibbo-Kobayashi-Maskawa (CKM) matrix Vqkm = VE^l^ ■ We 
furthermore introduce the notation 

f = —Y U Y U \ B = —Y D Y D \ L = —Y L Y L \ (18) 

In order to reconstruct the results for a general Yukawa sector, we multiplied each Feyn- 
man diagram by a factor (nh) m , where m denotes the number of fermion loops involving 
Yukawa couplings. After analyzing the structure of the diagrams that can arise, we could 
establish the following set of replacements that have to be performed in order to take into 
account a generalized Yukawa sector 

n h a t —> trT, n h ab —> trB, 

n h a T — > trL, nhO? t — > trT 2 , 

rihOtl — > trB 2 , rihCttCtb tiTB, 

rihCt 2 - — > trL 2 , n h a t ~^ (trT 1 ) 2 , 

n 2 h a 2 b -> (trB) 2 , n\a 2 T -> (trL) 2 , 

n1a t ab —> trTtrB, n 2 a t a T —> trTtrL, 
n 2 ai,a T — > trBtrL. (19) 



13 



Of course, only traces over products of Yukawa matrices can occur because they arise from 
closed fermion loops. Using Eqs. (IT7|) and (118]) it is straightforward to see that in Eq. (fl9|) 
only traces of diagonal matrices have to be taken except for trTB which is given by 



trTB = tr 




CKM 




(20) 



The addition of a fourth generation of fermions to the SM particle content can be also 
easily accounted for by this general notation. In this case, the Yukawa matrices become 
4x4 dimensional. If we assume that the fourth generation is just a repetition of the existing 
generation pattern but much heavier and if we neglect all SM Yukawa interactions, then the 
explicit form of Yukawa matrices reads 



3 x3 

«f 



with F 



T,B,L. 



(21) 



Here T and B stand for the up- and down-type heavy quarks, and L for the heavy charged 
leptons, while dp denotes the corresponding Yukawa couplings as defined in Eq. Q. Since 
in our calculation no Yukawa couplings for neutrinos have been introduced we cannot in- 
corporate heavy neutrinos. This would require a dedicated calculation which, however, does 
not pose any principle problem. 

We are now in the position to present the results for the beta functions of the gauge 
couplings which are given by 
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In the above formulas n G denotes the number of fermion generations. It is obtained by 
labeling the closed quark and lepton loops present in the diagrams. 

To obtain the results for the three-loop gauge beta functions, one also needs the one-loop 
beta functions of the Yukawa couplings, cf. Eq. ([6]). They can be found in the literature, 
of course. Nevertheless we decided to re-calculate them as an additional check of our setup. 
For completeness, we present the analytical two-loop expressions which read 
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The one-loop results have been expressed in terms of Yukawa matrices since these expres- 
sions enter the three-loop beta functions. At two-loop order, however, we refrain from 
reconstructing the general expression which would require an extension of the rules given in 
Eq. (US]). 

Our independent calculation of the two-loop Yukawa beta functions is also interesting as 
there is a discrepancy between [20J and [23] concerning the absence of terms proportional 
to a b a t X in Eqs. ( 12"5"1) and (12"6"1) . We were able to confirm the results in Ref. (23| . 

In Appendices \M and [B] we provide the results for the renormalization constants which 
lead to the beta functions discussed in this Section. 



IV. CHECKS 

We successfully performed a number of consistency checks and compared our results with 
those already available in the literature. We describe these checks in detail in this Section. 

The consistency checks show that all computed renormalization constants are local, that 
the renormalization constants of the gauge couplings are gauge parameter independent and 
that the beta functions are finite. We also find that the beta functions calculated by con- 
sidering different vertices in Lorenz gauge agree among themselves and with the results of 
the computation in BFG. 
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In order to test that the program FeynArtsToQ2E correctly translates FeynArts model 
files into model files for QGRAF/q2e, we reproduced the beta function for the Higgs self- 
coupling to one-loop order and the beta functions for the top and bottom quark, and the tau 
lepton Yukawa couplings to two-loop order (cf. previous Section). We have even considered 
quantities within the Minimal Supersymmetric Standard Model, like the relation between the 
squark masses within one generation, which is quite involved in case electroweak interactions 
are kept non-zero. Furthermore, we find that the divergent loop corrections to the BBB 
vertex vanish in the Lorenz gauge, as expected since for this vertex no renormalization is 
required. We performed the latter check up to three-loop order. 

Another check consists in verifying that in the vertex diagrams no infra-red divergences 
are introduced although one external momentum is set to zero. We do not have to consider 
two-point functions since they are infra-red safe. One can avoid infra-red divergences by 
introducing a common mass for the internal particles. Afterwards, the resulting integrals 
are evaluated in the limit q 2 3> m 2 where q is the non- vanishing external momentum of the 
vertex diagrams. This is conveniently done by applying the rules of asymptotic expansion [4l[ 
which are encoded in the program exp. The setup described in Section III CI is particularly 
useful for this test since exp takes over the task of generating FORM code for all relevant 
sub-diagrams which can be up to 35 for some of the diagrams; this makes the calculation 
significantly more complex. In our case the asymptotic expansion either leads to massless 
two-point functions or massive vacuum integrals or a combination of both. The former are 
computed with the help of the package MINCER, for the latter the package MATAD is used. As 
a result one obtains a series in m 2 /q 2 where the coefficients contain numbers and \n(m 2 /q 2 ) 
terms. For our purpose it is sufficient to restrict ourselves to the term (m 2 /q 2 ) and check 
that no logarithms appear in the final result. With this method we have explicitly checked 
that the W1W2W3 and three-gluon vertices are free from infra-red divergences. Since the 
results for the gauge coupling renormalization constants agree with the ones obtained from 
the other vertices also the latter are infra-red safe. 

Let us finally comment on the comparison of our findings with the literature. We have 
successfully compared our results for the two-loop gauge beta functions with [9| and the 



two- loop Yukawa beta functions with [23j. Even partial results for the three-loop gauge 
beta functions were available in the literature [l6| . That paper comprises all hitherto known 
three-loop corrections in a general quantum field theory based on a single gauge group, 
however, the presentation of the results relies on a quite intricate notation. Its specification 
to the SM is straightforward, however, a bit tedious. For convenience, the translation rules 



needed to specify the notation of [16| to ours is given in Appendix [Ql 

After specifying the notation of Ref. [lfj to ours we find complete agreement, taking into 
account the following modifications in Eq. (33) of 16[ :^ The terms 



7g 2 tT (Y a Y h ) tr (Y b Y a C(R)) 

g 2 tT (yqyfr) tr {y b Y c )C{S) ca 

12r ^ ' 

have to be "symmetrized" , so that they read 

l/ 7g 2 ti(Y a Y b )ti(Y b Y a C(R)) 7g 2 tr (Y a Y b ) tr (Y b Y a C(R)) 
2 I 127 127 



4 We want to thank the authors of [la ] for pointing out the issue of symmetrization and for assistance in 
deriving the translation rules. 
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1 ( 9 2 tr {Y a Y b ) tr (Y b Y c ) C(S) ca _ gHr (Y a Y b ) tr (Y b Y c ) C(S) Cu 
+ 2 I 12r 12r ' 

Furthermore, one has to correct the obvious misprint 

5c/ 4 tr (C(R) 2 Y a Y b ) 5# 4 tr (C{RfY a Y a ) 

+ 127 + 127 ' ^ 



V. NUMERICAL ANALYSIS 

In this Section we discuss the numerical effect of the new contributions to the gauge 
beta functions. We solve the corresponding renormalization group equations of the gauge 
couplings numerically and take into account the contributions from the Yukawa couplings 
and the Higgs self-coupling to two-loops. As boundary conditions we choose 

af^{M z ) = 0.0169225 ± 0.0000039, 
af 5 (M z ) = 0.033735 ± 0.000020 , 
af S {M z ) = 0.1173 ±0.00069, 
«P(M Z ) = 0.07514, 
ap(M z ) = 0.00002064, 
a^(Mz) = 8.077- 10" 6 , 

4ttA = 0.13, (31) 

where the first six entries correspond to experimentally determined values while the value 
for the Higgs coupling is determined assuming a Higgs boson with mass 125 GeV: 

4,aA 1252 G f .0.13. (32) 
2« 2 2-(V2-174) 2 GeV 2 

Note that the values in Eq. ( 13"T|) are given in the full SM, the top quark being not integrated 
out. For a description how these values are determined from the knowledge of their directly 



measured counterparts [40] , we refer to [40|, |42 



It is noteworthy that for all three gauge couplings the sum of all three-loop terms involving 
at least one of the couplings a&, a T or A leads to corrections which are less than 0.1% of the 
difference between the two- and three-loop prediction. 

In Fig. [3] the running of the couplings ax, «2 and a$, is shown from \i = Mz up to high 
energies. At this scale no difference between one, two and three loops is visible, all curves 
lie on top of each other. 

The differences between the loop orders can be seen in Fig. H] which magnifies the in- 
tersection point between a\ and ct 2 . There is a clear jump between the one- (dotted) and 
two-loop (dashed) prediction. The difference between two and three loops (solid curves) is 
significantly smaller which suggests that perturbation theory converges very well. 

The experimental uncertainties for a\{Mz) and ai{Mz) as given in Eq. fl3T|) are reflected 
by the bands around the three-loop results. Defining the difference between the two- and 
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FIG. 3. The running of the gauge couplings at three loops. The curve with the smallest initial 
value corresponds to a±, the middle curve to a 2 , and the curve with the highest initial value to 03. 



three-loop result as theoretical uncertainty one observes that it is smaller than the experi- 
mental one, however, of the same order of magnitude. Without the new three-loop calcu- 
lation performed in this paper the theory uncertainty is much larger than the experimental 
one. 

Also in the case of 0:3 perturbation theory seems to converge well. However, in contrast 
to <x\ and a 2 the experimental uncertainty turns out to be much larger than the theoretical 
uncertainty. This is not surprising as the relative experimental uncertainty of 0:3 at the 
electroweak scale is quite large compared to its electroweak counterparts. The relative 
experimental and theoretical uncertainty of is plotted as a function of the renormalization 
scale in Fig. [5j Note that by construction we have that Aa3/a3| t h e0 r y approaches zero for 
H ->• M z . 

Let us finally identify the numerically most important contributions. This is done by 
running from fj, = M z to n = 10 16 GeV and by comparing the contribution of an individual 
term to the total difference between the two- and three-loop prediction. Similar results are 
also obtained for lower scales. 

About 90% of the three-loop corrections to the running of the gauge couplings arises from 
only a few terms. In the case of a,\ there is only one term which dominates, namely the one 
of order a\a\. For a 2 one has a contribution of +56% from the O{a\o^) term, +13% from 
order a\otz and +37% from All other terms contribute at most 5% and partly also 

cancel each other. Except for the term of 0(af) all these terms are presented in this paper 
for the first time. 

The beta function /3s is dominated by the strong corrections, however, large cancellations 
between the a\ (+137%), a\a 2 (+45%), a\a 2 t (+28%) and a\a\ (+17%) terms on the one 
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FIG. 4. The running of the electroweak gauge couplings in the SM. The lines with positive slope 
correspond to a\, the lines with negative slope to a.2- The dotted, dashed and solid lines correspond 
to one-, two- and three-loop precision, respectively. The bands around the three- loop curves 
visualize the experimental uncertainty. 



hand and the a\a t (—112%) and a^a 2 a t (—16%) on the hand are observed. It is worth 
noting that the four-loop term of order a 5 s amounts to —58% of the difference between the 
two- and three-loop predictions. This number has been obtained by adding the four-loop 
QCD term [13, EH to 3 . 



VI. CONCLUSIONS 

In this paper we present three-loop results for renormalization constants that are used to 
compute the three SM gauge coupling beta functions, taking into account all contributions. 
We have checked that our expressions agree with all partial results present in the literature. 
Furthermore the two-loop corrections to the Yukawa couplings have been computed. We 
have performed the calculation using both Lorenz gauge within the unbroken phase of the 
SM and BFG in the broken phase. Our final result is valid for a generic flavour structure 
with an arbitrary CKM matrix. It is furthermore sufficiently general to consider a fourth 
generation of quarks and leptons. 

In order to perform the calculation in an automated way we have written an interface, 
FeynArtsToQ2E, between the package FeynArts and our chain of programs (QGRAF, q2e, 
exp, MATAD, MINCER) allowing to handle the O(10 6 ) Feynman diagrams, which have to be 
considered in the course of the calculation, in an effective way. Thus, we could perform 
several checks involving various different Green's functions. FeynArtsToQ2E is not limited 
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FIG. 5. Comparison of the relative experimental and theoretical uncertainty of 03. The theoretical 
uncertainty is given by the dashed line, the solid curve corresponds to the experimental one. 



to the SM but can easily be used for extensions like supersymmetric models. 



ACKNOWLEDGMENTS 

We thank David Kunz for tests concerning FeynArtsToQ2E and Mikolaj Misiak for sug- 
gesting to generalize the Yukawa structure of the final result. This work is supported by 
DFG through SFB/TR 9 "Computational Particle Physics". 



Appendix A: Renormalization constants 

In this Appendix we present analytical results for the renormalization constants Z ai , Z a 



a 3 i Zb-i Z\V, %Gi Zh, ^c w , ^c g i 



Z CG , Z cCW , Z cCG , Z www , Z GGG , Z HHW , where the definition 



of Z ai is given in Eq. (J3]) and the field renormalization is defined through 
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-W 
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Z B B 
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bare 



Z Cw C\V 



„bare 



Z W W 

Z CG C G 



G bare = Zq q ? 
^bare = ^jj 



(Al) 



B, W, G, cw and c G denote the gauge boson and ghost fields. The scalar field H is defined 
in Eq. f ]D5j) . The renormalization constants for the three-particle vertices are also defined 
in a multiplicative way. 

Some of the results listed below contain the gauge parameters £b, £w and They are 
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conveniently denned via the corresponding gauge boson propagator which is given by 
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with X = B, W, G. Note that £ x — 1 corresponds to Feynman and £ x = to Landau 
gauge. Our analytical results read 
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Note that in the results for Zh and Zhhw there are terms which contain explicitly a& 
and a t since we have not been able to reconstruct the corresponding expressions in terms of 
B and T. They drop out in the final result for Z a2 . 



Appendix B: Two-loop Yukawa coupling renormalization constants 

This Subsection contains the two-loop results for the Yukawa coupling renormalization 
constants defined through 
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Appendix C: Beta functions for oqed and sin 2 9yy 

This Appendix contains explicit results up to three-loop order for the QED coupling oqed 
and the weak mixing angle. We refrain from providing expressions for the renormalization 
constants but directly list the beta functions. They are obtained in a straightforward way 
from Eq. (J2]) and are given by 
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Appendix D: Comparison with Ref. [1 



In this Appendix we provide the explicit form of the expressions needed for the comparison 
with Ref. 16[. In this paper two-component Weyl spinors were used. To make contact with 



our convention based on four- component Dirac spinors we define 



D 



(Dl) 



where £ and \ are left-handed Weyl spinors and denotes a Dirac spinor. Thus, the 
Lagrange density of the SM can be expressed in terms of 45 Weyl spinors: 



Xt ) St ) Xb ) Q> ) Xt > £t j Xvt 5 Xc i £c j Xsi £s) Xfii Xvy.1 Xui 6uj Xd, £,d, Xe, Sej Xv a 



(D2) 
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For simplicity we have suppressed the SU (3) color indices for all quark spinors. Of course, 
each quark spinor has to be understood as a triplet in color space. In this basis the Yukawa 
matrices become 45 x 45 dimensional. 



In the notation of [16j] the part of the Lagrange density describing the Yukawa couplings 
is given by 



(D3) 



where Y a and Y a are the (complex conjugated) Yukawa matrices, <p a are real scalar fields, 
ip l and ip 1 are (Hermitian conjugated) spinor fields. There are four real scalar fields in the 
SM, which means that we have four Yukawa matrices. They are given by 
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Y 4 = -^ 



1 



U3x3 


U3x3 


Do o 

U3x3 


1 n /, Tl „ 


Do -, 
U3xl 


U3xl 


U3xl 


Do o 

u 3x3 


Do o 

u 3x3 


1 1 1 i Tl n n 


u 3x3 


U3xl 


Do 1 
U3xl 


Do ! 

U3xl 


03x3 


l?/t -11-3x3 


03x3 


03x3 


03x1 


03x1 


03x1 


h/fel-3x3 


3 x3 


3 x3 


3 x3 


3 xi 


3 xl 


3 xl 


0lx3 


0lx3 


0lx3 


0lx3 











0lx3 


0lx3 


0lx3 


0lx3 










0lx3 


0lx3 


0lx3 


0lx3 











V 



\ 



(D4) 



In the above formulas, all matrix elements not explicitly given are zero. Y 1 /Y 3 is the 
Yukawa matrix of the real/the imaginary part of the isospin down component of the SM 
Higgs doublet. Y 2 /Y 4 is the Yukawa matrix of the real/the imaginary part of the isospin 
up component of the SM Higgs doublet. So we have 
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2 + i0 4 



± 

v/2 V0 1 + i0 3 



(D5) 



The part of the Lagrange density describing the Higgs self-interaction reads in the nota- 
tion of Ref. [H 



(D6) 



In the SM we have 



Kaaa = 6 X (4ttA), 
Kabb = Kbab = Kbba = 2 X (47rA) (o / 6), 

Kbcd = (otherwise) . 



(D7) 



The expressions given above are generic for all three gauge groups. However, there are 
some expressions which depend on the specific gauge group one wants to consider. The re- 
mainder of this Section lists these expressions. For each of the three gauge groups, we give the 
expressions for the generators in the representations of the scalar fields, S A , and of the Weyl 
spinors, R A . We also give expressions for the invariants T(S),C(S),T(R),C(R),C(G),r, 
all of which are symbols used in Ref. [16]. They are defined as 



Tr (S A S B ) = 5 AB T(S), 
Tr (R A R B ) = 5 AB T(R), 
f ACD f BCD = 5 AB C(G), 



qA qA f~i( Q\ 

D ac°cb — ° \ D )ab, 



R Ak R Aj 

rAA 



C(R) 
r . 



(D8) 



j abc are ^ e structure constants of the respective gauge group. The indices A,B,C take 
values in the range 1, . . . , r, where r gives the dimension of the group. 
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1. [/(!) 



For the case of U(l) gauge group, we have 



R 1 = Diag 
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(D10) 



where the ellipsis is to be replaced twice by the first 15 entries. The derivation of S l is given 
below. The entries in R l are the hypercharges Y of the respective spinors. They can be 
derived by using the relation Y = Q — 1%, where Q corresponds to the electric charge and 
I3 is the weak isospin. Furthermore, we have 



T(5) = l, 

T(R) = 10 
C{G) = 0, 
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Let us now we show how to derive Eq. (ID9I) . We want to find the U(l) representation 
transforming the four real scalar fields of the Higgs doublet of the SM (see Eq. (ID5I) . The 
matrix S 1 is the generator of this transformation, 
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(1 + iuS 1 + O(oo 2 )) 



^3 



(D12) 



with u being the transformation parameter. In a next step we take advantage of the fact 
that it is known how the SM Higgs doublet transforms in order to determine S 1 . As the SM 
Higgs doublet has hyper charge 1/2, we have 



1 
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(D13) 



With the help of the last equation one can determine the transformation of the SM Higgs 
doublet, like, e.g., <p[ = (f) 1 — ^dj 3 + 0{uj 2 ). It is then straightforward to determine S 1 by 
inserting the equations found in this way in Eq. flD12j) . 
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2. SU(2) 



For the SU (2) group, we have 



/o -l\ 
0-10 
10 

\1 J 



s 2 = - 





/o 


-i 





o\ 


i 


i 











2 











-1 




V° 





1 


°) 



s 3 = -- 



and 



1/2 





3 X3 


° 1^2^3x3 


3 X3 


3 xl 


3 xl 


3 xl 


3 x3 


3 x3 


3 x3 


3 x3 


3 xl 


3 xl 


3 xl 




3 x3 




3 x3 


3 xl 


3 xl 


3 xl 


3 x3 


3 x3 


3 x3 


o 3x3 


3 xl 


3 xl 


3 xl 


0lx3 


0lx 3 


0lx3 


0lx3 


°2,2 





°2,1 


0lx3 


0lx3 


0lx3 


0lx3 











0lx3 


0lx3 


0lx3 


0lx3 


°1,2 





<1 



( 1 \ 
-1 
-10 
\o 10 J 

\ 



\ 



(D14) 



(D15) 



In the last equation a A are the Pauli matrices. Not all of the matrix elements left out in 
R A vanish, but these elements play no role for the comparison of our results to [l6j]. The 
derivation of the generators S A proceeds in analogy to the derivation of S 1 explained in 
the former Subsection. One merely has to substitute the generator of £7(1), |l, by the 
generators of SU(2), \a A , in Eq. (IDT3]) . 
We furthermore have 



T{S) = 1, 

T(R) = 6, 
C(G) = 2, 



C(S) — -l4 X 4, 

C(R) = Diag (1, 1,1,0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 1, 
r = 3 . 



(D16) 



The ellipsis has to be replaced by the first entries twice. 



3. SU(3) 

In the SM the matrices S A vanish for the group SU(3). The matrices R A are block- 
diagonal and read 



R A = - BlockDiag (A A , -(A A ) T , A A , -(A A ) T , 0, 0, 0, ... ) 



(D17) 



The ellipsis has to be replaced twice by the former entries, which contain the Gell-Mann 
matrices A A . 

Finally, we have 



T(S) = 0, 



C(S) — 04x4, 
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T(R) = 6, C(R) = - BlockDiag (l 12xl2 , 3x3 , Ii2xi2, 3x3 , 1 12x12, 3x3 ) , 

C(G)=3, r = 8. (D18) 
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